0 


* 


AD-A204  378 


REPORT  DOCUMENTATION  PAGE 


U  REPORT  SECUROY  CLASSIFICATION 

UNCLASSIFIED 

Vo  HESTkw:Vi>rt  MAkifAihs  ^ 

2«  SECURITY  CLASSIFICATION  AUTHORITY 

3  OISTRIBUTION/AVAILABIUTY  OF  REPORT 

26  OECLASSIFiCAriON  DOWNGRADING  SCHEDULE 

Approved  for  pubuc  release;  distribution  is  unlimited. 

4  PERFORMING  ORGANIZATION  REPORT  NUMBERlSt 

5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6«  NAME  OF  PERFORMING  OROANtZATlON 

00  OFFICE  SYMBOL 

74  NAME  OF  MONITORING  ORGANIZATION 

tit  bpoftcMtht 

Naval  Ocean  System*  Center 

NOSC 

Naval  Ocean  Systems  Center 

6c  ADDRESS  I'C'ff'  St*to  9na  ZiP  Codei 

7b  ADDRESS  f’C'fV.  Swe  ondZJP  Code) 

San  Diego,  CA  921S2-S000 

San  Diego,  CA 

92152-5000 

04  NAME  Of  FUNDING  SPONSORING  ORGANIZATION 

80  OFFICE  SYMBOL 

9  PROCUREMENT  INSTRUMENT  lOENTlFlCAHON  NUMBER 

lit  tCfihcMilOl 

Space  and  Naval  Warfare 

SPAWAR 

Systems  Command 

8c  ADDRESS -Oiy.  Srbfo  ood  ZfP  Codo; 

10  SOURCE  OF  FUNDING  NUMBERS  | 

program  element  NO 

PROJECT  NO 

TASK  NO 

AGENCY 

REWSON  System*  Program 

ACCESSION  NO 

Washington,  DC  20363-S100 

RDDA 

MP57 

NSA 

DN288  591 

1 1  TITLE  Sfcufity  CI»Ssihe9tioni 

Accuracy  of  MUF  and  Field  Strength  Prediction*  by  HFBC84 


12  personal  autmop-s* 


R.A.  SprAgue 


type  9f  report 

Fro(««iioaAl  Paper 


13b  TIME  COVERED 

May  1987  May  1987 

conu  *  TO  ' 


14  DATE  OF  REPORT  yttr  Month  D»ri 

November  1987 


’6  SUPPLEMENTARY  notation 


1  1  7  COSATi  COOES  | 

FIELD 

GROUP 

SUB  GROUP 

50  Subject  TERMS<Coni>r»t/ea"'?v«r$«»Toec«#Ar’y^^/^?rif>fy&y6foc*ouml 


frequency  range 
field  itrength 
median  ralue* 


a 


5  9  ABSTRACT  Contmug  on  re^orsm  ‘f  notosson/  Offti  itjontrfy  Ov  biock  ni/mo*r/ 


DTIC^ 

^fE6 


A  Btatletical  analysii  of  the  accuracy  of  MUF  and  field  etrength  prediction  by  HFBC84,  a  prediction  program  developed 
for  the  HF  Broadcast  World  AdministraCive  Conference,  Is  presented.  A  database  of  13,054  hours  of  oblique  sounder  MOFs 
measured  on  70  paths  was  used  to  obtain  the  residuals  in  the  predicted  MUFs.  The  MUF  model  had  a  bias  of  1.17  MHe 
low,  an  rms  error  of  4.67  MHt,  and  a  correlation  coefficient  with  the  measured  data  of  0.83.  In  the  4000  to  5000  km 
range,  the  MUF  model  has  a  large  negative  bias  (predicts  high).  It  has  Us  highest  bias  during  summer  when  it  is  2.2  MHs 
low.  During  winter  it  has  bias  of  1.1  MHs  low.  The  model  reaches  a  minimum  in  bias  at  1800  local  time  at  the  path 
midpoint;  whereas,  Us  maximum  error  Is  offset  12  hours  from  the  minimum.  A  modified  version  of  CCIR  Database  C  was 
used^  to  obtain  residuals  in  predicted  field  strength;  only  the  8l  paths  tor  short  path  propagation  were  retained  (12,277 
median  values).  The  results  of  the  study  are  presented  as  a  function  of  the  ratio  of  frequency  to  predicted  MUF.  In  the 
f/MUF  range  from  0.4  to  1.0  the  bias  in  predicted  field  strength  goes  from  2.4  dB  low  to  3.1  dB  high  with  a  peak  value 
at  3.8  dB  low  at  f/MUF  -0.7.  In  this  same  frequency  range,  the  rms  error  goes  from  12.7  dB  to  13.2  dB.  In  the  f/MUF 
range  firom  1.0  to  1.5,  the  bias  goes  from  3.1  dB  high  to  1.5  dB  low,  and  the  rms  error  goes  from  13.2  dB  to  13.0  dB, 
peaking  at  17.0  dB  at  f/MUF  -1.3.  Beyond  this  range,  both  the  bias  and  rms  error  increase,  rapidly  reaching  a  bias  of 
39.6  dB  low  and  an  rms  error  of  40.5  dB  at  f/MUF  -2.0.  These  results  indicate  that  the  field  strength  method  used  in 
HFBC84  for  over-the-MUF  frequencies  should  be  re-examined. 
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A  scaciscical  analysis  of  Che  accuracy  of  NUF  and  field  screngch  prediccion  by  HF3C34,  a 
prediccion  program  developed  for  the  HF  Broadcasc  World  Adminiscraciva  Conference,  is 
presenced.  A  daca  base  of  13,054  hours  of  oblique  sounder  ttOFs  measured  on  70  pachs  was  used  zo 
obtain  Che  residuals  in  the  pradicced  MUFs.  The  MUF  model  had  a  bias  of  1.17  MHz  low,  an  rms 
error  of  4.67  MHz,  and  a  correlation  coefficient  with  Che  measured  daca  of  0.33.  In  the  4000- 
5000  km  range,  the  MUF  modal  has  a  large  negative  bias  (predicts  high).  It  has  its  highest 
bias  during  sunmer  when  it  is  2.2  MHz  low.  During  winter  it  has  a  bias  of  1.1  MHz  low.  The 
model  reaches  a  jinimum  in  bias  ac  1800  local  cime  at  Che  path  mid-poinc;  whereas,  its  maximu.T. 
error  is  offset  12  hours  from  Che  minimum.  A  modified  version  of  CCIR  Daca  Base  C  was  used  to 
obtain  che  residuals  in  predlcced  field  strength;  only  Che  81  pachs  for  shore  path  propagation 
were  retained  (12,277  median  values).  The  results  of  che  study  are  presented  as  a  function  oi 
che  ratio  of  frequency  co  predicted  MUF.  ..In  Che  f/MUF  range  from  0.4  to  1.0  the  bias  in 
predicted  field  screngch  goes  from  2.4  dB  low  cm,  3.1  dB  high  with  a  peak  value  ac  3.3  dB  Low  a: 
f/MUF  -  0.7.  In  this  same  frequency  range,  the  rlss  error  goes  from  12.7  dB  to  13.2  dB.  In  the 
f/MUF  range  from  1.0  to  1.5,  che  bias  goes  from  3,1. dB  high  co  1.5  dB  low,  and  the  rns  error 
goes  from  13.2  dB  Co  13.0  dB,  peaking  ac  17.0  dB  ac  f/MUF  -  1.3.  Beyond  this  range,  both  the 
bias  and  rms  error  increase,  rapidly  reaching  a  bias  of  39.6  dB  low  and  an  rms  error  of  40.5  dB 
ac  f/MUF  -  2.0.  These  results  indicate  chat  che  field  strength  method  used  in  HFBC84  for  over- 
the-MUF  frequencies  should  bo  re-examined. 


The  effective  operation  of  long  distance  high  frequency  (HF)  communications  has  increased 
in  proportion  co  the  ability  to  predict  variations  in  the  Ionosphere.  These  variations  are  af¬ 
fected  in  a  complex  manner  by  solar  activity,  seasonal  and  diurnal  changes,  as  well  as  latitude 
and  longitude.  Such  a  predictive  capability  has  permitted  communicators  Co  optimize  fre¬ 
quencies,  antennas  and  ocher  circuit  parameters. 

Initially,  manu;  1  methods  were  developed  for  analyzing  ionospheric  variations  on  HF  cir¬ 
cuits  of  short,  intermediate,  and  long  distances.  Because  Che  manual  methods  were  laborious 
and  Cime  consuming,  various  organizations  have  developed  computer  programs  co  analyze  HF  cir¬ 
cuit  performance. 


One  such  prediccion  method  is  Chat  adopted  by  che  First  Session  of  che  World 
.■administrative  Radio  Conference  (WARC)  for  HF  Broadcasting  (1984)  .  This  computer  program  is 
commonly  called  HFBC84  .  Simultaneously,  Recommendation  COM  5/1  (1984)  was  adopted  which  re¬ 
quested  che  CCIR  Co  undertake  some  further  studies  including  studies  concerning  this  method 
In  particular,  it  requested  che  CCIR  co  provide  daca  necessary  to  refine  ~wo  numerical  con¬ 
stants  used  in  the  field  screngch  method.  The  approach  taken  by  the  CCIR  was  to  (1;  have 
Interim  Working  Party  6/1  collate  measurement  data  of  hourly  median  field  strengths  available 
during  1984  and  (2)  establish  a  new  Interim  Working  Patty,  6/13,  co  undertake  the  task  o: 
determining  the  numerical  constants.  The  field  screngch  measurement  daca  collated  by  IVP  6.1 
is  known  as  Daca  Base  C.  This  daca  was  used  by  IWP  6/13  to  determine  che  two  numerical  con¬ 
stants  in  the  field  screngch  method;  these  constants  were  subsequenclv  adopted  bv  the  CCIR  i;-. 
1986  , 


This  paper  presents  che  uncertainty  assessment  of  che  predicted  maximum  usable  frequency 
(MUF'i  and  field  screngch  from  HFBCSa.  The  version  tested  contains  the  new  numerical  constants 
as  recommended  by  che  CCIR.  The  method  of  HF  model  uncertainty  assessment  starts  with  the  con¬ 
struction  of  che  daca  bases  percinonc  to  the  parameter  being  assessed.  For  each  model  the 
observed  and  modeled  daca  ara  compared  using  a  daca  screening  program  (Sailors  ec  al..  1981. 
Sailors  ec  al.,  1936).  Error  cables  are  produced  as  a  function  of  various  relevant  propagation 
parameters . 
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DATA  BASES  USED  IN  THE  ANALYSIS 


The  observed  dace  base  used  in  Che  HUE  analysis  eons  is  Cs  of  13,054  median  hourly  hOF 
values  derived  from  70  dlfferenc  oblique  sounder  paths.  Path  lengths  range  from  196  km  to 
10.376  km  vith  a  relatively  uniform  distribution  among  intermediate  lengths.  The  data  Is  also 
approximately  uniformly  distributed  in  season,  local  time  and  12  month  running  mean  sunspot 
number  (1959-1981).  It  contains  data  from  nearly  all  geomagnetic  north  latitudes  vith 

relatively^Iess  at  low  latitudes  and  very  little  at  southern  latitudes. 

The  field  strength  data  base  is  a  modified  version  of  Data  Base  C.  Our  version  of  Data 
Base  C  consists  of  12,277  median  skywave  field  strength  data  values  (dB  above  IfiV/m)  normalized 
to  1  kU  erp  measured  on  31  short  paths.  This  differs  from  Data  Base  C  in  that  all  "long  path" 
data  has  been  removed. 

The  data  making  up  our  version  of  Data  Base  C  was  collected  worldwide  over  the  years  1964- 
1984  and  contains  R12  values  in  the  range  13-163.  In  most  cases  several  frequencies  were  used 
on  each  circuit  which  allows  investigation  of  the  frequency  dependent  aspects  of  the  model. 
Unfortunately,  no  frequencies  below  4.8  MHz  are  represented  so  the  field  strength  variation 
with  frequency  near  the  LUF  could  not  be  adequately  tested. 

ANALYSIS  PROCEDURE 

The  data  screening  computer  program  used  to  analyze  the  residuals  between  the  predicted 
and  obser'zed  data  is  called  0ASCR3.  DASCR3  is  able  to  partition  a  large  data  base,  at  runtime, 
into  sub-data  bases  each  of  which  contains  only  chat  data  which  satisfies  a  preset  screening 
criteria  (e.g.  ,  path  length  between  1000  and  2000  km).  Because  this  partitioning  is  done  at 
runtime,  it  is  not  necessary  to  physically  retain  multiple  sub-data  bases. 

The  analysis  procedure  consists  of  running  the  respective  model  programs  to  generate  a 
predicted  data  base  which  matches  all  Che  conditions  of  Che  observed  data.  Then  DASCS3  is  used 
to  determine  the  statistics  of  Che  residuals  for  Che  sub-data  base.  Statistics  are  determined 
for:  (1)  residual  -  observed  datum  -  predicted  value;  (2)  relative  residual  - 

residual/observed  datum;  and  (3)  absolute  relative  residual  -  absolute  residual/obser'ved  dac-jm. 

KUF  CALCULATION  IN  HFBC84 

The  F2  layer  basic  HUF  calculation  in  HF5C84  is  based  on  numerical  maps  of  the  F2  layer 
ordinary  wave  critical  frequency,  foF2,  and  the  obliquity  factor  for  a  3000  km  path,  M(3006), 
For  paths  less  chan  4000  km  the  F2(4000)I1UF  is  determined  as  the  product  of  fot2  and  .MCIOOO)  at 
Che  path  midpoint  and  the  numerical  factor  1.1,  interpolating  in  R.j,  as  needed  The  F2  layer 
MUF  for  these  paths  is  Chen  found  by  interpolation  between  TY (0.)  MUF  (i.e.,  foF2)  and 
F2(4000)MUF  and  subsequent  multiplication  by  an  M-facCor  which  consists  of  a  polynominal  in 
path  length. 

For  paths  greater  Chan  4000  km  a  simple  control  point  method  utilizing  Che  path  midpoint 
and  control  points  2000  km  from  each  terminal  is  used.  The  F2(4000)MUF  is  determined  at  each 
of  these  control  points  and  the  minimum  of  these  values  is  taken  as  the  path  F2  layer  .MUF 

For  paths  less  than  4000  km  Che  basic  median  E- layer  .MUF  is  also  calculated.  In  this  case 
an  empirically  derived  calculation  of  foE  given  by  Davies  (1965). 

foE  -  0.9  ((180+1.44  R.j^2>  cosx'l®'^^  (MHz), 

is  used  instead  of  numerical  maps.  Here  x'  is  a  simple  function  of  the  calculation  point 
zenith  angle  (Davies,  1963).  A  control  point  method  similar  to  the  above  is  used  in  this  case 
•■•ith  points  1000  km  from  each  terminal  used  for  paths  greater  chan  2000  km.  The  path  MUF  is 
determined  as  the  maximum  of  the  values  determined  for  Che  F2  and  E  layers. 


MUF  ANALYSIS  RESULTS 

Results  of  our  comparison  for  the  entire  data  base  are  summarized  in  Table  I.  Some  o: 
these  results  can  be  misleading  in  that  they  may  hide  a  large  deviation  for  a  particular  sec  o: 
parameters  which  is  compensated  by  a  similarly  large  variation  in  the  opposite  sense  for 
another  sec.  Overall  however,  performance  of  the  MUF  prediction  scheme  in  HF3C34  is  quite 
good,  with  an  average  residual  of  1.17  MHz,  a  rms  deviation  of  4,52  .Mhz  and  an  overall  correla¬ 
tion  coefficient  of  .827.  Results  for  the  ocher  statistics  are  also  shown  in  Table  1. 

As  shown  in  Table  1,  Che  MUF  prediction  capability  of  HFBC34  is  quite  good.  Results  for 
two  problem  areas  for  which  additional  modeling  effort  may  be  necessary  are  presented  below. 


Figure  1  shows  the  average  residual  as  a  function  of  sounder  path  length.  Also  indicated 
is  Che  standard  deviation  for  each  case.  This  'igure  indicates  Chat  HFBC84  has  a  marked  ten¬ 
dency  to  over  predict  Che  KUF  for  paths  in  the  4000-5000  km  range  and  under  predict  in  the 
7000-8000  kn  range. 

Inspection  of  Che  data  base  shows  chat  a  significant  majority  of  the  data  in  the  4000  km- 
5000  km  range  is  less  chan  4500  km  and  most  is  less  Chan  4300  km.  This  is  approximately  the 
transition  distance  for  a  2-hop  mode  change.  Similarly,  the  7000-8000  km  range  is  just  less 
than  Che  distance  for  Che  next  mode  change.  These  results  taken  together  indicate  that  the 
method  used  for  prediction  of  Che  KUF  at  these  critical  distances  is  not  adequate. 

For  paths  S  4000  tos,  Che  three  control  points  are  approximately  coincident.  Choice  of  the 
minimum  as  Che  path  MUF  Ignores  Che  effect  of  the  mode  switch  in  reducing  the  MUF. 
Determination  of  a  hop  length  dependent  M-factor  in  these  cases  should  produce  the  required 
reduction  in  Che  tflIF  at  these  ranges. 

Similarly,  for  paths  £  8000  km  Che  H-factor  has  become  quite  large.  Ignoring  its  effect 
in  the  MUF  determination  causes  the  under  prediction  seen  in  Figure  1.  Proper  inclusion  of  an 
M- factor  at  distances  greater  than  4000  tos  should  produce  improved  predictions  in  this  range. 

In  Figure  2  we  show  the  average  residual  as  a  function  of  month  (season)  .  Again  we  also 
indicate  the  standard  deviation  for  each  case.  This  figure  indicates  chat  the  model  has  a  ten¬ 
dency  Co  under  predict  in  summer  and  winter  months.  The  correlation  coefficient,  shown  in 
Figure  3,  also  shows  a  reduction  in  summer  months. 

Figure  4  shows  the  average  residual  and  the  corresponding  standard  deviation  as  a  function 
of  midpach  local  tine.  The  minimum  error  occurs  for  the  hours  1500  to  1900  hours  LT.  The  max¬ 
imum  bias  occurs  from  0500  to  1100  hours  LT,  an  average  offset  of  12  hours  from  the  fflinimu.m 
error.  RFBCS4  has  its  lowest  standard  deviation  from  1800  to  0700  hours  LT. 

The  seasonal  and  diurnal  errors  in  HF3C84  probably  relate  to  the  polynominal  expression 
used  for  the  M  factor.  The  expression  is  only  a  function  of  path  length  and  hence  must  assume 
a  constant  layer  height  independent  of  season  and  time  of  day  in  the  derivation. 


Table  1.  70  path  statistical  analysis  summary  for  .HFBC84  MUF  prediction. 


STATISTIC 


ANALYSIS  RESULT 


total  path  hours  13054 

average  residual  1.17  MHz 

rms  residual  4.67  MHz 

average  relative  residual  .059 

rms  relative  residual  .  242 

average  absolute  relative  residual  .217 

correlation  coefficient  .827 


Figure  I.  Average  residual  (bias>  and  standard  deviation  of  the  predicted  KU?  as  a  function  of 
path  length. 
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HELD  STRENGTH  CALCUIATION  IN  HFBC84 

The  field  strength  calculation  in  HFECS^  is  divided  into  path  length  regimes.  One  method 
is  used  for  path  lengths  less  than  7000  kn;  another  method  is  used  for  path  lengths  greater 
than  9000  km;  and  from  7000  km  Co  9000  km  Che  field  strength  is  determined  by  interpolating  as 
a  function  of  range  the  values  of  the  field  strengths  at  7000  km  and  9000  km. 

For  paths  less  Chan  or  equal  to  7000  km,  Che  calculated  field  strength  depends  on  whether 
the  operational  frequency  is  in  relationship  to  the  path  MUF.  For  frequencies  less  than  or 
equal  to  Che  MUF,  the  field  strength  calculation  Includes  terms  for  basic  free  space  loss, 
ionospheric  absorption,  ground  loss,  auroral  absorption,  transmitter  power  and  antenna  gain, 
and  a  numerical  factor  7.3  dB  included  to  account  for  those  effects  of  skywave  propagation  not 
otherwise  included.  The  ionospheric  absorption  is  Che  same  as  used  by  the  CCIR  (1970).  The 
nighttime  ionospheric  absorption  is  due  to  the  work  of  Vakal  (1961).  A  loss  of  2  dB  is  as¬ 
signed  for  each  ground  reflection.  A  cable  of  experimentally  determined  factors  is  used  to 
represent  auroral  absorption.  The  Cable  depends  on  path  length,  time,  season,  and  geomagnetic 
latitude.  The  values  peak  at  polar  latitudes  with  values  from  3  colO  dB.  For  frequencies 
greater  chan  the  MUF,  an  additional  loss  mechanism  known  as  over-che-HUF  loss  is  added.  this 
is  included  Co  account  for  Che  fact  chat  observed  field  strength  does  not  abruptly  vanish  as 
the  operating  frequency  begins  to  exceed  the  basic  MUF  for  a  path.  In  HFBC85  the  loss  term  is 
given  by  130  [)f/MUF)  and  is  not  allowed  to  become  larger  chan  85  dB. 

For  ranges  greater  chan  or  equal  to  9000  km,  the  HFBC84  field  strength  method  uses  the  FTZ 
method  (Beckman,  1967;  DamboldC,  1975).  Ac  these  path  lengths,  the  effects  of  ionospheric 
tilts  and  scattering  from  irregularities  cause  contributions  to  Che  measured  field  strength 
from  paths  outside  Che  main  great  circle  path.  Beckman  determined  that  the  field  strength 
rises  from  a  frequency  (f^)  for  which  the  field  strength  is  some  low  limiting  value  to  a  peak 
at  some  intermediate  frequency  somewhat  below  Che  MUF,  and  then  decreases  to  some  low  limiting 
value  at  a  frequency  greater  chan  the  MUF  which  he  calls  the  operational  MUF  (fo) .  The  f,  -fo 
combination  then  determines  the  transmission  range  of  frequencies  for  which  usable  field 
strength  may  exist.  The  lower  frMuency.  fj^.  not  a  LUF,  is  the  frequency  for  which  the  field 
strength  is  0  dB  (*iV/m)  f^r  10  Tctf  effective  radiated  power  or  the  frequency  at  which  free- 
space  field  strength  for  10'’  kW  effective  radiated  power  equals  the  ionospheric  absorption 
assuming  no  ocher  loss  terms.  The  parameter  f.  is  a  function  of  a  season,  the  12  month  running 
mean  sunspot  number,  angle  of  incidence  at  90  iSm  assuming  mirror  reflection  from  300  km  and  hop 
span  0  for  great  circle  path  ray,  solar  zenith  angle  at  a  height  of  90  km  number  of  legs  (twice 
the  number  of  hops),  winter  anomaly  term,  and  slant  path  length  d.  For  F2-layer  modes  fo  - 
K’MUF  where  K  is  an  empirically  derived  factor  and  is  a  function  of  the  MUF  for  Che  hour  of  in¬ 
terest,  the  maximum  and  minimum  MUF  for  Che  2A-hour  day,  and  a  sec  of  constants  dependent  on 
raypach  azimuth.  The  constants  used  in  HFBC8^  have  been  changed  from  earlier  values  used  by 
FTZ,  and  continue  to  change  as  comparison  with  data  at  FTZ  warrants.  Given  then  f^^,  f<5.  f,  the 
transmit  power  and  antenna  gain,  the  field  strength  can  be  determined  from  a  simple  expression. 


FIELD  STRENGTH  ANALYSIS  RESULTS 

Results  of  our  comparison  for  the  entire  data  base  are  summarized  in  Table  2.  The  results 
are  given  for  all  observation  frequencies,  for  observation  frequencies  less  than  equal  the 
predicted  .MUF,  and  for  observation  frequencies  above  the  MUF.  Clearly,  there  is  a  considerable 
difference  in  accuracy  of  HFBC84  field  strengths  between  frequencies  less  chan  equal  the  MUF 
and  frequencies  above  Che  MUF.  It  is  because  of  this  large  difference  that  the  rest  of  the 
discussion  of  the  field  strength  prediction  accuracy  is  divided  according  to  observation  fre¬ 
quencies  relative  to  the  MUF. 

Again  as  was  Che  case  in  the  MUF  accuracy  results  only  the  problem  areas  can  be  high¬ 
lighted.  Figure  5  shows  Che  avet  -ge  residual  and  its  standard  deviation  as  a  function  of  the 
obser'.'ation  frequency  f  Co  .MUF  ratio.  Below  the  .MUF,  the  model  is  quite  accurate:  yet.  it 
predicts  somewhat  high  as  the  frequency  approaches  the  MUF.  It  continues  to  predict  high  until 
a  ratio  of  about  1.5  is  reached.  At  ratios  of  15  or  larger  the  model  predictions  are  clearl;.- 
not  satisfactory. 


Figure  6  presents  the  average  residual  as  a  function  of  circuit  length  for  both  fre¬ 
quencies  less  than  equal  the  MUF  and  greater  than  the  MUF.  For  frequencies  less  than  equal  the 
:-iuF,  HFBC84  predicts  high  for  short  paths  (0-1000  km),  predicts  moderately  high  for  paths  from 
2000  km  CO  7000  km,  and  beyond  7000  km  performs  poorly.  For  frequencies  greater  than  the  .MUF. 
HFBC84  predicts  low  for  path  ranges  from  0  to  2000  km,  high  for  path  ranges  2000  km  to  7000  km. 
and  low  for  path  lengths  greater  than  7000  km. 
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Figure  7  shovs  Che  seesonel  resulce  for  frequencies  less  chan  equal  Che  HUF.  HFSC84  shows 
a  large  bias  occurring  during  sumser.  Ac  frequencies  above  che  HUF,  HFBC84  Is  approximacely 
unlforo  over  all  seasons  In  nose  eases. 

Al chough  Che  sunspoc  behavior  foe  HFBC84  ac  frequencies  less  chan  equal  che  HUF  shows  an 
Increase  In  bias  with  Increasing  soooched  sunspoc  number,  Che  effecc  of  varying  sunspoc  numbers 
Is  noc  as  large  as  Che  Increase  In  bias  for  frequencies  greacer  chan  che  HUF  as  figure  8  shows . 

Finally,  figure  9  shows  che  resulcs  as  a  funecion  of  midpach  local-clme.  Koce  chac  for 
frequencies  less  Chan  equal  Co  Che  HUF,  che  residual  Is  uniform  chrough  che  day;  whereas,  for 
frequencies  greacer  chan  che  HUF,  che  bias  Is  a  minimum  ac  local  noon  with  very  high  bias  ac 
nighc. 


Table  2.  81  pach  scaclsclcal  analyses  summary  for  HFBC84  field  screngch  predlccion. 


STATISTIC 

cocal  pach  hours 
average  residual 
rms  residual 

average  relaclve  residual 
rms  relaclve  residual 
average  absoluce  relaclve  residual 
correlaclon  coefflclenc 


OVERALL 

f<HUF 

f>HUF 

12277 

8783 

3494 

3.97  dB 

1.93  dB 

9.10 

dB 

16.68  dB 

12.06  dB 

24.73 

d£ 

1.32 

0.04 

4.53 

29.59 

18.22 

47.35 

3.21 

2.34 

47.35 

.75 

.78 

,57 

DISCUSSION  AND  RECOHMENDATIONS 

The  previous  secClons  presenced  che  resulcs  of  a  comparison  of  che  accuracy  of  che  .HUF  and 
field  screngch  predicced  from  HFBC84.  FuncCional  discrepancies  were  noced  in  chose  seccions. 
In  chis  section  possible  causes  for  che  discrepancies  are  discussed  and  recommendations  for  im- 
provemencs  are  made.  The  HUF  and  field  screngch  models  are  created  separately. 

For  che  HUF  computation,  two  problem  areas  were  noced.  The  first  dealt  with  the  fact  chat 
che  .HUF  calculation  Ignores  Che  effect  of  minimum  hop  mode  changes  with  increasing  range.  The 
correction  is  simple;  Insert  the  range  of  one  hop  of  a  multihop  mode  into  che  H- factor  equa¬ 
tion.  This  should  reduce  che  bias  occurring  ac  ranges  where  mode  changes  occur.  The  seasonal 
and  diurral  errors  probably  relate  Co  che  fact  chac  the  polynomlnal  expression  used  for  the  .M- 
faccor  Is  only  a  funcclon  of  paCh  length  and  hence  muse  assume  a  constant  layer  height  and 
semi- thickness  independent  of  season  and  time  of  day.  Lockwood  (1983)  has  presenced  a  non- 
iceracive  procedure  which  enables  evaluaclon  of  Che  HUF  using  a  distance  factor,  with  allowance 
for  variations  in  both  peak  height  and  changes  in  underlying  plasma.  If  the  expression  for 
mirror  height  of  reflection  given  by  Lockwood  (1984)  for  frequencies  at  the  HUF  is  utilized,  it 
would  be  possible  to  determine  more  precisely  the  effects  of  change  in  number  of  hops  on  the 
MUF . 


In  interpretating  che  field  screngch  results,  it  should  be  realized  that  unlike  ionos¬ 
pheric  characteristics  such  as  che  HUF  or  foF2  which  are  scaled  from  ionograms  in  che  same  way 
throughout  che  world,  procedures  used  to  obtain  field  strengths  vary  in  both  method  and  in 
quality  Usually  the  input  voltage  at  Che  receiver  input  in  mv  is  measured.  This  measured 
value  is  Chen  converted  to  field  strength  based  on  some  knowledge  of  losses  inherent  in  both 
receiving  and  transmitting  systems.  In  addition,  certain  assumptions  ere  made  for  a  nominal 
takeoff  angle  for  the  received  signal.  When  a  data  base  such  as  Data  Base  C  is  assembled, 
there  is  some  inherent  inaccuracies  in  it  which  cannot  be  assessed. 


Aside  from  these  considerations,  the  overall  performance  of  che  field  strength  model  is 
quite  poor  and  could  be  improved,  particularly  in  its  over-che-HUF  prediction  method.  It  is 
recommended  that  the  full  Wheeler  (1966)  model  be  implemented  in  HFBC84  for  ranges  loss  than 
"000  km.  This  would  give  the  needed  range,  diurnal,  and  sunspoc  variations  noc  included  in  the 
present  method.  For  short  paths  for  frequencies  less  chan  equal  che  HUF,  HFBCSi*  predicts  high 
This  might  be  due  in  part  to  the  method  used  to  determine  takeoff  angle.  HFBC34  uses  an  e;;- 
pression  from  Lockwood  (1984)  which  Lockwood  states  is  valid  for  frequencies  0.7  to  0.95  times 
the  .HUF  and  not  valid  for  short  paths.  The  problem  with  che  use  of  che  expression  employed  at 
short  range  is  chat  che  mirror  height  of  reflection  is  no  longer  linear  with  decreasing  range, 
Ac  short  range  the  mirror  height  of  reflection  is  similar  to  chat  for  frequencies  near  the  HUF. 
onlv  noc  as  pronounced.  Sailors  et  al .  (1986)  provides  a  correction  to  this  expression  for 
ranges  less  than  D„..^  which  is  also  given  therein.  Winter  anomaly  corrections  to  ionospheric 
absorption  such  as  cMc  due  to  Schultz  and  Gallet  (1970)  might  improve  its  seasonal  accuracy. 
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Figure  8.  Average  residual  (bias)  of  Che  predicted  field  screngch  as  a  function  of  smoothed 
sunspot  number. 
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Figure  9.  Average  residual  (bias)  of  the  predicted  field  strength  as  a  function  of  raid-path 
local  time. 

For  ranges  greater  chan  7000  km,  %>here  the  FTZ  method  is  employed,  Sprague  (1987)  found 
chat  Advanced  PROPHET  was  more  accurate  than  HFBC8A,  particularly  at  frequencies  greater  than 
MUF.  Both  programs  employ  the  FTZ  method.  However,  Advanced  PROPHET  does  not  use  the  FTZ  ex¬ 
pression  for  operational  MUF.  It  uses  1 . 85(HPF)MUF ,  where  HPF  is  the  90  percentile  value 
assuming  a  normal  MUF  distribution.  For  an  HPF  of  1.15  this  is  the  98.3  percentile  point.  At 
f/MUF  greater  than  2.0,  Advanced  PROPHET'S  bias  also  rapidly  increases  reaching  a  peak  at  f/MUF 
-3.3.  This  might  imply  that  2.87(HPF)MUF  be  used  for  the  operational  MUF  which  would  be 
equivalent  to  the  99.95  percentile  point  of  a  normal  distribution.  The  trouble  with  the  ex¬ 
pression  used  in  the  FTZ  method  for  operational  MUF  is  chat  it  is  based  on  paths  terminating  in 
Germany.  The  result  is  Chat  it  weights  the  European  area.  In  addition,  FTZ  changes  Che  con¬ 
stants  as  measurements  indicate.  A  method  such  as  employed  in  Advanced  PROPHET  with  perhaps  a 
higher  percentile  point  chan  98.3  has  more  universal  application. 
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